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ABSTRACT 

We present post-outburst observations of the mid-infrared spectrum and sub- 
millimeter continuum of the illuminating source of the newly-discovered McNeil's 
Nebula in the L1630 region of Orion. The 12 /mi flux of this source has increased 



by a factor of ~25 after the outburst, whereas the submillimeter continuum 
00 . 

remains at its pre-outburst level. The bolometric luminosity has increased by 

O. 

■^j- , mid-infrared spectrum exhibits a strong and red continuum with no emission or 



absorption features. The infrared slope of the spectral energy distribution char- 
acterizes the illuminating source as a flat-spectrum protostar, in both its active 
and quiescent states. New CO spectral line observations show no evidence for a 



O, 

molecular outflow. 

CO 



Subject headings: stars: formation — stars: pre-main sequence — stars: variables 



1. Introduction 

Recently, McNeil (2004) has reported the discovery of a bright reflection nebula in a 
previously dark region of the L1630 cloud within the Orion B complex. An infrared source at 
the apex of McNeil's Nebula (hereinafter referred to as McNeil's Nebular Object, or MNO) 
has significantly brightened in the optical and near-infrared since 2003 November. Detailed 
imaging and spectroscopy of that object are described by Reipurth & Aspin (2004), Briceno 
et al. (2004), and Vacca et al. (2004). Apparently the nebula has brightened in the past, 
indicating that MNO has resumed its outburst activity after roughly 37 years of quiescence 
(see the discussions by Reipurth & Aspin 2004; Briceno et al. 2004). 

In the low state, MNO is positionally coincident with a faint /-band point source 
(Eisloffel & Mundt 1997), a 2MASS source (J05461313-0006058), an IRAS point source 
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(05436—0007), an ISO source (Abraham et al. 2004), as well as a clump of submillimeter 
emission known as LMZ 12 (Lis, Menten, & Zykla 1999) and OriBsmm55 (Mitchell et al. 
2001). Previous work has shown that the L1630 cloud contains a number of deeply embedded 
young stellar objects (YSOs), which drive powerful outflows (see Simon et al. 2004, and ref- 
erences therein). In the immediate vicinity of MNO there are several Class (HH24-MMS, 
HH25-MMS) and Class I (SSV59, SSV63, HH26-IR) sources, as well as a series of outflows 
(HH 22-27). Both before (Eisloffel & Mundt 1997; Lis et al. 1999) and after (Reipurth & 
Aspin 2004) the current episode of activity, it has been suggested that MNO is the origin of 
the HH23 jet. In this paper, we provide new mid-infrared and submillimeter observations of 
MNO in its high state and we discuss their implications for understanding the evolutionary 
status of that object. 



2. Observations and Data Reduction 

Submillimeter continuum photometry at 450 and 850 /im of MNO was obtained on UT 
2004 March 10 with the Submillimeter Common User Bolometer Array (SCUBA) on the 
James Clerk Maxwell Telescope (JCMT) at Mauna Kea, Hawaii. The precipitable water 
vapor (PPV) level was steady at ~1.8mm, with frequent skydip observations showing zenith 
optical depths of 0.36 at 850 /im and 2.05 at 450 /mi. The data were acquired in 20 sets of 
18 s integrations in a small 9-point jiggle pattern with the secondary mirror chopping 60" in 
azimuth at 7.8 Hz. Mars and CRL618 were observed to check the telescope pointing and to 
determine the absolute flux calibration. The pointing has an rms error of <2"; the flux cali- 
bration is accurate to ~10% at 850 //m and ~30% at 450 /mi. Standard reduction techniques 
were performed using the SURF software provided by the Joint Astronomy Center. 

Mid-infrared spectroscopic observations of MNO in the iV-band were carried out on UT 
2004 March 11 at the United Kingdom Infrared Telescope (UKIRT) on Mauna Kea with 
the Michelle imaging spectrometer (Glasse, Atad-Ettedgui, & Harris 1997). The Michelle 
detector is a Si: As 320x240 pixel array operating between 8 and 25 /mi. In spectroscopy 
mode, the instrument uses a 91" slit (in our case with a P. A. = 0° orientation). The Low-N 
grating and a 2-pixel wide (0 / .'76) slit were used to provide spectral coverage between 7.5 
and 12.5 //m at a resolving power of R ~ 200. Atmospheric variations were removed using a 
standard nod-and-chop pattern during the total on-source integration time of ~1250 s. The 
atmosphere was dry (PPV ~1.6mm) and stable, providing FWHM values for point sources 
of 0"8 at 10 /im. Slit images demonstrate that MNO is point-like in the mid-infrared. 

The spectra were sky-subtracted, flat-fielded and combined using the Observatory Re- 
duction and Acquisition Control (ORAC) software (Allan et al. 2002). The IRAF task 
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APALL was then used to extract spectra with a 3 pixel (1'.'14) aperture. A wavelength 
calibration solution was determined from terrestrial emission lines identified with ATRANS 
atmospheric models (Lord 1992). The rms of the wavelength fit was ~0.04yum. The primary 
flux standard Rigel (BS 1713) was used for both telluric correction and flux calibration. The 
MNO spectrum was divided by the standard star spectrum and the result was scaled by a 
model thermal spectrum based on the iV-band (10 fxm) flux of the standard star. 

Spectral line maps of 12 CO (J=2 — 1) and 13 CO (J=2 — 1) around MNO (covering 5'x5' 
and 2'x2', respectively) were obtained from Mauna Kea with the 230-GHz heterodyne re- 
ceiver on the CSO telescope on UT 2004 March 17. The observations were conducted in the 
on-the-fly mode, providing Nyquist-sampled maps with a spatial resolution of ~31". The 
50-MHz AOS backend provided a spectral resolution of ^O.OGkms^ 1 per channel and rms 
noise levels of ~1.1 K in both lines. Standard procedures with the CLASS software package 
were used to reduce the data. 



3. Results 

The broadband spectral energy distribution (SED) of MNO is shown in Figure 1 for both 
the high and low activity states. In the active state, the /-band flux is from Briceno et al. 
(2004, UT 2004 February 26), the additional optical (g'r'i') and near-infrared (JHK f ) data 
are from Reipurth & Aspin (2004), and the L and M-band fluxes are from the spectrum of 
Vacca et al. (2004). In the quiescent state, the data are the same as those used by Abraham 
et al. (2004) to analyze the pre-outburst SED, with the inclusion of an J-band flux point from 
Briceno et al. (2004, UT 1999 January 9). The standard photometric zero-points of Bessell 
(1979), Fukugita et al. (1996), and Tokunaga (2000) were used to convert magnitudes into 
flux densities. The error bars in Figure 1 represent the systematic uncertainties in absolute 
flux calibration. 

Consistent with the factor of ~16 near-infrared brightening of MNO, the 12 ^m flux 
density rose by a factor of ~25 during the high state. In the standard classification scheme 
for young stellar objects, or YSOs (Lada & Wilking 1984; Adams, Lada, & Shu 1987, Andre 
et al. 2000), a SED type is assigned from the logarithmic slope of the spectrum, n, where 
vF v oc v n . For the conventional infrared indices of the post-outburst SED of MNO, we 
obtain n = —0.62 ± 0.14 between 2 and 10 /im (corresponding to the K — N color) and 
n = —0.66 ± 0.13 between 2 and 12 jum 1 . The 2MASS 2yum and IRAS 12 /im points from 



lr These indices are calculated from flux measurements taken ^5 weeks apart, and are therefore subject to 
any changes in the near-infrared fluxes during that time. 



-4- 



the quiescent phase of MNO yield an index of n — —0.35 ± 0.18. The foregoing values of the 
SED slope in the infrared establish MNO as a flat-spectrum/Class I object regardless of its 
activity level. 

Figure 2 shows the mid-infrared spectrum of MNO from 7.5 to 12.5 //m. The spectrum 
exhibits a bright, red, featureless continuum across the iV-band. A simple power-law fit 
to the continuous spectrum gives F v oc A 2 ' 5 , which corresponds to n = —1.5 ± 0.1. The 
positions of various solid-state features that are commonly seen in the spectra of other YSOs 
are labeled in Figure 2, although none of them are detected here. 

The flux densities at 450 and 850 fim measured in our SCUBA observations of MNO are 
1.589±0.099 and 0.316±0.005 Jybeam -1 (quoted errors are la statistical uncertainties), and 
correspond to FWHM apertures of ~9" and ~14", respectively. Our 850 /im flux density is 
essentially identical to the value cited by Mitchell et al. (2001) for a ~20" diameter aperture, 
based on SCUBA maps of L1630 made by those authors in 1998. Lis et al. (1999) have 
argued that the 1.3 mm emission is confined within a deconvolved diameter of < 3". Given 
the agreement between the 850 /xm flux densities measured in two different sized apertures, 
the circumstellar structure responsible for the submillimeter emission must also be compact, 
and < 10 4 AU in diameter if we assume the accepted distance of d ~ 450 pc for L1630. 

Contrary to the behavior of the optical and infrared spectrum of MNO, no changes are 
apparent in the submillimeter continuum (350 yum < A < 1.3 mm) as a result of the current 
activity. Assuming the Rayleigh- Jeans limit and optically-thin emission, the submillimeter 
flux density follows the relationship F v oc z/ 2+/3 , where (3 is the power-law index of the grain 
opacity. A least-squares fit to the four submillimeter points in the (log z/, log F v ) plane gives 
(3 = 0.65±0.19. The slightly more sophisticated model of Mitchell et al. (2001, their Equation 
2), if applied to the 450 and 850 /im measurements, yields a range of (5 between 0.79 and 
1.24 (for dust temperature values T dust ~ 50 and 20 K, respectively). The discrepancy in [3 
between a simple fit and the Mitchell et al. (2001) model may indicate that the material is 
not optically-thin and/or the Rayleigh- Jeans approximation is inappropriate. The value of 
j3 is apparently near unity, although we caution that it would be best determined with a full 
SED model fit if more far-infrared data become available while MNO is in an active state. 



4. Discussion 

As shown above and discussed by Abraham et al. (2004), the small (slightly negative) 
infrared spectral index, n, establishes MNO as a flat-spectrum/Class I object. The large 
spectral slope in the submillimeter is also more typical of Class I sources than the presumably 
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more evolved Class II sources (Dent, Matthews, & Ward-Thompson 1998). It is worthwhile 
to compare MNO with the other young outflow sources in L1630, since it is reasonable to 
assume these other objects, as siblings of MNO, may experience similar outbursts. Using 
our data and the flux values given by Mitchell et al. (2001), we note that the 850 /im flux 
decreases by a factor of 4 between the Class (e.g., HH 24-MMS and HH 25-MMS) and 
Class I (e.g., SSV 59 and SSV 63EW) objects in L1630, and by another factor of 4 from the 
Class I objects to MNO. With an appropriate correction for distance, flat-spectrum sources 
in other clouds have 850 ^m fluxes similar to that of MNO (e.g., HL Tau). Applying the 
standard assumption of optically-thin, isothermal dust («850 — 0.03 cm 2 g _1 , T^ ust = 50 K; 
Beckwith 1999), we infer from the 850 jum flux a total (gas + dust) circumstellar mass of 
M cs ~ 0.06 M . Our mass estimate is an order of magnitude lower than that calculated 
by Abraham et al. (2004), primarily due to the different T^ ust assumed. Nevertheless, M cs 
is significantly larger than the values typically found for Class II disks. The submillimeter 
evidence therefore agrees well with the interpretation from the infrared SED that MNO 
is currently in transition from the Class I to Class II stage of protostellar evolution. The 
circumstellar environment of MNO is probably dominated by a massive accretion disk, with 
only a remnant envelope structure. 

The mid-infrared spectrum of MNO shown in Figure 2 is notably devoid of any solid- 
state spectral features. None of the three major PAH features (at 7.7, 8.6, and 11.2 /im) are 
seen, in agreement with the absence of the more frequently detected 3.3 xim PAH feature in 
the 1-5 /im spectrum presented by Vacca et al. (2004). PAH features are seldom observed in 
low-mass YSOs (Smith, Sellgren, & Tokunaga 1989; Brooke, Sellgren, & Geballe 1999), but 
are more common in their higher mass counterparts, the HAeBe stars (Ressler & Barsony 
2003; Habart, Natta, & Kriigel 2004). There is no evidence for 3.3 xim PAH emission in the 
infrared spectra of the other Class I objects in L1630 discussed by Simon et al. (2004). More 
puzzling is the absence of a 9.7 ^m amorphous silicate feature in the spectrum of MNO. 
According to Cohen et al. (1984), at least two of the other Class I objects in L1630 (SSV 59 
and 63) have strong silicate absorption bands. The optical depth of the 3.1 /im H2O ice band 
measured by Vacca et al. (2004), r ice ~ 0.7, suggests that the 9.7/im feature of MNO should 
also be in absorption. Using information provided by Brooke et al. (1999, their Figure 8), we 
would expect a band depth of r(9.7/im) ~ 0.4. Such a shallow absorption could be entirely 
filled in by emission from the optically-thin surface layers of a disk if it is viewed at a suitably 
high inclination angle (Chiang & Goldreich 1999, their Figure 5). 

Recent VLA observations at wavelengths of 3.6, 6, and 20 cm detected no radio emission 
from MNO to rms noise levels of 42, 48, and 140yuJy, respectively (Claussen 2004). The flux 
densities expected from dust, based on an extrapolation from the submillimeter wavelengths, 
fall below these detection limits. Of the Class I sources in L1630, radio detections have been 
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reported previously for SSV63 (three separate components: Reipurth et al. 2004) and HH26- 
IR, whereas SSV59 has a 3a upper limit of 66/iJy (Gibb 1999). The radio emission from 
objects like SSV59 and MNO is either absent or optically-thick. 

The lack of an obvious high velocity jet or molecular outflow from MNO remains to be 
explained and further distinguishes this object from the other YSOs in L1630. Although a 
powerful wind is evident in post-outburst optical and infrared spectra, shocked emission from 
the standard optical and infrared forbidden lines is absent (Reipurth & Aspin 2004; Briceho 
et al. 2004; Vacca et al. 2004). The submillimeter CO spectral line maps obtained by us and 
by Lis et al. (1999) show no spatial distinction between redshifted and blueshifted emission. 
The limited spatial resolution of these maps does not definitively rule out the presence of a 
molecular outflow, and so interferometric observations will be required to address this issue 
appropriately. However, as shown in Figure 3, our single-dish observations do rule out the 
presence of gas moving at velocities greater than ~3km s" 1 from the rest velocity of the 
L1630 cloud. The best candidate outflow signatures from MNO are the HH23 clumps of 
[S II] emission, which are located ~3' to the north, along the direction of McNeil's Nebula 
(Eisloffel & Mundt 1997; Lis et al. 1999; Reipurth & Aspin 2004). 

Using a simple trapezoidal integration of the SEDs shown in Figure 1, we estimate that 
the bolometric luminosity of MNO has changed from ~3.5 to 34 L Q during the outburst. The 
pre-outburst value of Lboi was estimated by Lis et al. (1999) to be 2.7 L Q and by Abraham 
et al. (2004) to be 5.6 L . The former result is identical to ours when allowance is made for 
the different distances that were used. The Abraham et al. (2004) value is slightly larger 
than ours because those authors chose to correct for an assumed (non-local) extinction of 
Ay = 13 mags. However, such a correction leads to an over-estimate of L^i if the extinction 
is due to local circumstellar material (as is thought to be the case for MNO) because the 
short-wavelength flux is thermalized and then re-emitted at longer wavelengths. Our post- 
outburst value of L^i likely underestimates the true value because we lack information in 
the SED near its peak in the far-infrared. We have used a linear extrapolation between the 
infrared and submillimeter portions of the SED to estimate a peak flux of ~80 Jy (at roughly 
70 /im), which leads to an upper limit on the post-outburst L bo i < 90 L . 

Briceho et al. (2004) have calculated an intrinsic luminosity of the post-outburst MNO 
of L = 219 L Q from their /-band data, assuming that the spectral type is A0V and that 
Ai = 7.2, or Ay ~ 15. However, the post-outburst colors of MNO suggest a lower extinction 
of Ay ~ 11 (Reipurth & Aspin 2004), which is in good agreement with the value indicated by 
the depth of the infrared H 2 ice band (Vacca et al. 2004). Using the interstellar extinction 
law of Mathis (1990), the main sequence colors and bolometric corrections of Kenyon & 
Hartmann (1995), and the lower post-outburst extinction value, the Briceho et al. (2004) 
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luminosity is reduced by more than a factor of 4, to L^ ~ 47 L Q (this assumes Ry = 
3.1; the value would be roughly a factor of 2 larger if Ry = 5.0). The revised value is 
thus in good agreement with the estimate obtained by integrating under the SED. Both 
analyses, all the same, should be treated with some caution. In the post-outburst stage, 
there will be a considerable contamination of the /-band magnitudes from scattering in 
McNeil's Nebula (particularly with the large pixel scale in those observations). Moreover, 
it is not at all clear that main sequence colors and bolometric corrections apply to YSOs 
like MNO. At the same time, a simple integration of the SED includes any luminosity 
from accretion, and consequently may not be an accurate representation of the intrinsic 
photospheric flux (including the portion that is thermalized by local dust and re-emitted at 
longer wavelengths). With those caveats in mind, the L^oi values both before and after the 
outburst do not support earlier claims that MNO is a B type star, to order of magnitude 
they are similar to those of other YSOs in L1630 (Cohen et al. 1984; Berrilli et al. 1989), 
and they are significantly lower than those derived for FU Orionis stars (e.g., Sandell & 
Weintraub 2001, their Table 1). 



5. Summary 

We have presented new observations in the mid-infrared and submillimeter of the out- 
burst star which illuminates McNeil's Nebula. The object has brightened by a factor of 
~25 in the mid-infrared, yet remained at the same brightness in the submillimeter. The 
bolometric luminosity has increased by (at least) an order of magnitude, but still remains 
low compared to that expected from either an early-type photosphere or a FUor. Our lim- 
ited CO spectral-line maps indicate that no high-velocity flows of molecular gas have yet 
appeared in the immediate vicinity of this object. Far-infrared observations of MNO from 
the Spitzer Space Telescope would help to complete the post-outburst SED and aid in deter- 
mining Lb i and Td ust , thus M cs , providing important evolutionary constraints on this object. 
Observations of solid-state features in the mid- and far-infrared spectrum of MNO would 
also serve to provide a constraint on the disk inclination angle and thus would offer a test 
of our explanation for the absence of the 9.7 /zm feature. 

We are grateful to Tom Kerr, Sandy Leggett, Thomas Lowe, Sandrine Bottinelli, and 
Jonathan Williams for their assistance with the observations and data reduction. Suggestions 
from an anonymous referee have greatly improved this Letter. We also want to thank Mark 
Claussen, Mike Cushing, Joel Kastner, Bo Reipurth, and Bill Vacca, who have generously 
shared their results before publication. 
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Fig. 1. — The broadband SED of MNO. Open circles mark data taken in the quiescent state 
before 2003 November, and filled circles mark post-outburst data from 2004 February and 
March (see the text for references). The mid-infrared spectrum is shown in greater detail in 
Figure 2. The continuum slope in the infrared indicates that MNO is a flat-spectrum/ Class 
I protostar. 
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Fig. 2. — The mid- infrared spectrum of MNO. A 5-pixel gap near 9.4 /im has been blanked 
out due to a residual telluric ozone line. Although no spectral features have been detected, 
the positions of some of the more common YSO features are indicated. A simple power-law 
fit to the continuum (F u oc A 2 ' 5 ) is shown for comparison. 
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Fig. 3. — Composite CO (J=2— 1) spectra from the region surrounding MNO. Each spectrum 
is the average in a 1' diameter aperture, labeled according to its offset in declination in 
arcminutes from the position of MNO. The spectra are offset in T& by 3K for clarity. The 
dotted lines mark velocity offsets of ±3 km s _1 from the rest velocity of the L1630 cloud. 



